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Abstract. Supercooled liquid is described in terms of the theory of defected states 
of bond orientational order using the continuous theory of defects. It is shown 
that the initiation of a local orientation order leads to the breakdown of the initial 
SO (3) symmetry of the theory, as a result the plastic deformation fields of the 
linear defects become Abelian. The expressions for linear dislocation and disclination 
tensor potentials are derived within the framework of the quasistationary adiabatic 
approximation. It is shown that a-relaxation and /3-relaxation of supercooled liquid 
can be described as relaxation processes in the disclination and dislocation subsystems 
accordingly. Disclination topological moment is proposed to use as an order parameter 
at the description of glass transition. 
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1. Introduction 

It is known that the translation symmetry absence in liquids (does not entail) (is 
not accompanied by) an absolutely random atom order. Most liquids under low 
temperatures possess a short-range or middle-range order. This ordering is determined 
by the directions and lengths of valence bonds, valence angles, extents of the atom 
and ionic radiuses, etc. In contemporary journals this order is called "local order" [1], 
"locally favored structures" [2] or "bond orientational order" [3j, while the liquid state 
is called Frenkcel phase [3], high density phase [5] and so on. The pioneering liquid state 
theory, based on the bond orientational order conception, was suggested by Frenkcel |6J. 
In spite of the objective complications of experimental investigations and theoretical 
description of liquids there is appreciable progress in development of this conception. 
The plenty experimental results, which sustain this conception, (see [4] for example) 
and theoretical works, which develop one, [7] have recently appeared. Among of them 
let us note the model of defected states of bond orientational order, suggested by Nelson 
[3] with the purpose of a disordered molecular systems description. Our interest to this 
approach is caused by the opportunity to use it for description of the supercooled liquids 
and the glass transitions [8]. 

The general idea of the considered model can be formulated in the following way: 
Under the high temperatures a liquid is an isotropic (quasi-gas) state of the condensed 
matter, which is characterized by absence of any short-range atom ordering and, as 
a result, locally SO(3)-symmetric. Under the temperature decrease or the pressure 
increase the initial system ordering happens. In this case the short-range order, which 
is locally symmetric relative to some point group, forms. Since this symmetry depends 
only on atom interaction, it is natural to assume that in most cases they are not among 
the Fedoroff groups, therefore, the ground state of these systems is a noncrystalline 
(a long-range order is absent) and infinite degenerated one owing to the geometric 
frustrations. 

The topological defects (dislocations and disclinations), which correspond to these 
frustrations, are basic structure elements of the system. It is important to emphasize 
that in contrast to the crystal defect, which increase the system energy, the defects in the 
considered model are essential structure elements, which compact the pacing and make 
the disordered system more energy-optimal than the crystal one [3]. In order to destroy 
the topology defects great energy consumptions are necessary. However, because of the 
frustrations and the thermal flotations the system of topology defects is moving. This 
motion is a cooperative atom motion, and due to this the condensed matter maintains 
a fluidity. 

2. The gauge field theory Lagrangian 

The gauge defect theory is a method of the disordered molecular system description, 
which is based on the mathematics of the gauge field physics. This approach implies a 
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fundamental role of the defects in the properties formation of the disordered molecular 
systems. Its basis was developed in the works of Lihkachev [9], Kadich and Edilen [10J, 
Nelson [3], etc. 

In the most fundamental form, which was suggested in [10], the gauge defect theory 
is represented within the generalized theory of the Yang-Mills gauge fields. To write 
down the Lagrangian of the defect system, let us make use of the standard theory of 
elasticity. The simplest Lagrangian, describing the system with elastic deformations, 
has the following form: 

1 i 

L o = ttAAxAXi - o ( M ^) 2 + 2 P u ij u ij] , 
z o 

where Xi{?i t) is the elastic strain field, A and \i are the Lame constants (the elasticity 
and shear modulus accordingly), p is the mass density (which is considered to be 
constant for simplicity), and the relative deformation components: 

Uij = Cij - 8ij = diXidjXi - <\j 

(Roman letters a, b, . . . are used to denote the space components set {i} = {1, 2, 3}, 
and Greek letters a, ft, . . . are used to denote the full index set, including the time 
component X 4 , {a} = {1, 2, 3, 4}). 

According to the gauge theory of dislocations and disclinations, the plastic 
deformation of the matter structure can be considered as a breakdown of the 
homogeneity of the rotation and translation (SO(3)>T(3)) groups action. In order 
to take into account these homogeneity breakdowns, the compensating fields (A% and 
(fl) are introduced into the Lagrangian, and the transition from ordinary to covariant 
derivatives is effected: 

^atXi ^ B a i 9 a Xi ~\~ £ijl-A a jXl (paii 

where are three generating matrixes of the semisimple group SO (3). After that, the 
Lagrangian L is replaced by the new Lagrangian 

L = Lq + s^Ly + saLa, 

where s v and sa are free parameters of the theory, and the first term describes elastic 
properties of matter: 

Lo = ^PoBuBi^ — - [X(Eii) 2 + 2fiEijEij] , 

where = BuB^ — 8ij is the unit strain tensor. The second term, 

1 

describes the dislocations; the following notation is used here: 

Di a /3 = d a <fi/3 — dpifia + Euj (Ai a ipj/3 — Aip(fj a + Ei a/3 Xj) , 



Fia/3 — d a Aip — dpA ia + EijiAjaAip, 



A description of a- and ^-relaxation supercooled in liquids 



4 



The third term, 

SaLa = —^A^ijlFijl + ^jSAFi^FiU, (1) 

describes disclintions. The Yang-Mills fields, A ia and <fi a , describe disclinations and 
dislocations accordingly. 

The gauge defects theory is a general one, since it considers the case of the ideal 
isotropic matter, which is locally symmetrical with respect to the SO(3)-group of 
rotations in the three-dimensional space. This theory is non-Abelian one (7J [10], which 
is the main problem of the considered approach. However, this paper will argue the idea 
that in the systems, which have a bond orientational order, this trouble disappears. 



3. The Lagrangian symmetry breakdown 

Note, that a molecular system can be SO(3)-symmetrical just as Hamiltonian only at 
the relatively high temperatures. At low temperatures, according to the considered 
model, a local bond orientational order appears there. Defects can exist only in local 
ordered systems. Then the system's choice of a ground state and appearance of an 
order parameter, Q(r), with some internal state space, M n , lead to the spontaneous 
symmetry breakdown. As a result the ground state symmetry does not agree with the 
Hamiltonian's SO(3)-symmetry. 

The broken symmetry is the bond orientational order symmetry, which is generally 
determined by the nature of the chemical molecular bonds. Let us consider two 
important cases: 1) In the systems with the spherical interaction potential the 
tetrahedral pacing is the most preferred one, which depends the icosahedral form of 
the short range ordering. It is known that the icosahedral symmetry group, Y, is 
isomorphous to the finite group of the proper rotations of M 3 space, which has three 
orbits (see appendix 1). Each orbit correspond to the rotation a round one of the 
symmetry axes. Since each of these orbits coincides with the corresponding internal 
state space, one can assert that the symmetry groups of these internal states are 
point symmetry subgroups in the SO (3). 2) Director is the order parameter in the 
nematic liquid. The ground state of this system can be arbitrary configuration of the 
codirectional directors. The system's «choice» of one of the ground states leads to the 
SO (3)— >SO(2) symmetry breakdown, as a result the system remains symmetrical in 
relation to rotations around the selected axis. 

Thus in the polytetrahedral and monoaxial nematic models the symmetry groups, 
which keep the order parameter, are subgroups of the SO (2). The nonlinear over Aij 
part of disclination field tensor vanishes, and the fields of the linear topological defects 
become Abelian: 

Fia/3 = d a Ai/3 — dpA ia + EijiAj a Aij3 — > F ial 3 = d a Ai/3 — dpA ia . 

This leads to a significant simplification of the considered systems analysis, since the 
Lagrangian is out of nonlinear A^ terms, which are typical for the Yang-Meals theory. 
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4. The linear topology defects fields in the quasistationary adiabatic 
approximation 

In order to find the tensor potentials of the linear defects one needs to do some 
simplifications: First of all let us describe the defects system in the quasistationary 
approximation (d^A ia = 0). In case of slow relaxation this seems to be appropriable. 
Besides, in this paper we will be interested in the typical for amorphous systems low- 
angled disclinations with v ~ 0.2 Frank index. In this case the considered model can be 
linearized: 

Xi = $iaXa + euf } + e 2 uf ] + . . . , 

A a i — zA a i + s 2 A^- + . . . , (2) 

9 (1) 

fod = etfai + £ Pod + • • • , 

where e ~ v is the small parameter, u t = euf 3 + e 2 uf S> + ... is the displacement vector. 
For simplicity let us assume that the free topology defects do not interact with each 
other. Below we will show that all these simplifications are equivalent to the adiabatic 
approximation. This seems to be appropriable in view of the significant difference 
between the dislocation and cisclination subsystem relaxation times. 

In order to determine the tensor potential of the linear disclination segment let us 
minimize the disclination term of the Lagrangian, 

In [H] it was shown that the quasistationary solution for the gauge field tensor potential 
of the linear disclination segment, dl a , has the form of 

A iA = x a ^ + 2C^j dl h A ik = EikjXj + C^j dk. (3) 

It is significant that in this expression there is no a summation over a, and C\ = 
according to A^(oo) = 0. 

Then from (fl]) one can determine the linear dislocation field. In the linear 
approximation the Atp product is proportional to e 2 , therefore one can neglect the first 
two terms in the brackets: 

where 

Fia/3 = d a Ai/3 — dpA ia . 

As it was noted above the linear approximation leads to the neglect of the dislocation 
birth-annihilation processes on the disclinations. 

The dislocation and disclination systems have intricate configurations. Therefore, 
in spite of the above approximations, the problem of the theoretical description of such 
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systems is a difficult task. But it can be significantly simplified when the high defect 
mobility is considered. In case of high defect concentration they screen the elastic 
energy of each other by the formation of the disclination and dislocation loops. The 
screening leads to the minimization of summary local topology charge. This minimizes 
the dislocation and disclination subsystems energy, since the elastic energy density of the 
defect loop decreases faster than one of the single linear defect. In case of disclinations 
it is 1/r 3 against 1/r. As a result the tensor potential of the disclination subsystem field 

is determined by the topology disclination moment, Qij(r): = Of (see appendix 
2). Substituting this potential form to the above expression one can write 

D ia/ 3 = d a (fip - dpifia + J e tt j {d a Ai P (x - f) - dpAi a (x - f )) x J dr = d a ip il3 - dpip ia . 

It follows that the linear dislocation field potential has the same form like the disclination 
one: 



Then one can get the tensor potentials of arbitrary configurations of the dislocation 
and disclination fields. Taking into account the screening the interaction energy of the 
linear defects system can be represented in the form of the interaction energy of the 
local topology moment system, and the problem is reduced to the known one. 

5. Topology moment interaction 

In spite of the fact that the linear dislocation and disclination tensor potentials have the 
similar form, they have different contributions in the elastic interaction. This determines 
the significant difference between the disclination and dislocation subsystem relaxation 
dynamics. 

The stress field is produced by the dislocations and disclinations. Keeping only 
nondiagonal elements of the deformation tensor in the first infinitesimal order one can 
write it in the form 

a ai ~ fiEai = n (BajBji - 5 ai ) = fi [e aa iXiA^ + e aa xiA" + (p ai + ip ia + diu a + d a Ui] . 
The elastic interaction energy, 



includes the inter-dislocation, U T , and inter-disclination, U R , interaction terms, as well 
as the dislocation-dislocation interaction U RT term. In spite of the fact that the last 
term is important in the calculations, in order to simplify the qualitative analysis we 
will limit ourselves only with the first two terms. 

As it was noted above, the elastic interaction energy of the disclination network can 
be represented in the form of the interaction energy of the topology moments system: 
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where a and (3 label the moments, a is the typical linear dimension of disclination loops, 
and 

The elastic energy of the dislocation topological moments system has the form of 



where 




Comparison of the expressions for the disclination interaction and dislocation interaction 
makes it clear that the first one is the long-range interaction 

J R (r) ~ 1 /r 3 , 

while the second one is short-range: 

J T (r) ~ 1 /r 5 , 5 > d = 3. 

This is a very important difference, since it leads to the difference of relaxation dynamics 
of the disclination and dislocation subsystems, which is the cause of the experimentally 
observed decoupling between rotational and translational diffusion [12] . 

In order to display it let us represent the energy functional of the disclination 
subsystem in the form of the functional of the frustration limited domains theory [71 [13]: 

U R = J {\\d»Q\x)\ 2 + ±T\Q R (x)\ 2 }d 3 x-± J J Q R (x)J R (\x-y\)Q R (y)d 3 xd\ 

where Q is some local structural variable (order parameter). The first item is a short- 
range interaction caused by continuity of Q, the second item is a weak long-range 
interaction, which describes the disclination moments interaction. In [T5] it was shown, 
that presence of a long-range interaction in such system leads to its frustration. It 
allows us to assert, that because of the long-range interaction the disclination subsystem 
is geometrically frustrated. The system's frustration leads to the slow non-Arrenius 
dynamics of relaxation processes proceeding in this system [7J, [14] . Thus in the 
considered system the relaxation of the disclination subsystem is a-relaxation. 

On the other hand the short-range interaction of the dislocation subsystem does 
not lead to the same frustration, and relaxation dynamics of this subsystem has a usual 
Arrenius form. It is /3-relaxation. 

This difference of the subsystems relaxation dynamics agrees with experimental 
results [12], which prove that in supercooled liquid there are two types of relaxation 
processes. These types concern the rotary and translational types of diffusion. The 
degrees of freedom, which correspond to disclinations moves, are inherently rotational 
and correspond to cooperative molecular motion, whereas the dislocation degrees of 
freedom are translational ones. 
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6. Topology moment as an order parameter 

In some theories [El [[3] the liquid-glass transition is considered as a special form of 
phase transition in the glasses and structurally disordered systems. In this case the 
problem is reduced to description of the interactive disclination subsystem. As in any 
another phase transition theory the central problem of this theory is the definition of 
order parameter. On the one hand an order parameter should characterize the local 
state and describe the short-range order symmetry, on the other hand the free energy 
of the system should be invariant under scaling transformations, so that the scaling 
hypothesis was true. 

The linear defect topology moment, Q^, used above satisfies these criterions. First 
of all this order parameter characterizes the short-range order symmetry by definition. 
At the same time it has the non-molecular (supermoleqular) scale, and the theory 
is scale-invariant. One can explain this in the following way: the local ordering of 
the topological moments under decreasing of the temperature leads to the symmetry 
breakdown, SO (3)— >SO(2). In this system the next generation of topology moments 
appears. At freezing in this system, just as in nematic, an ordering of the moments 
occurs. However, as a result of their random situation a general direction for all 
topological moments of the system does not exist. As a result this leads to the 
appearance of next generation of topological defects with a bigger scale. After changing 
the scale one can describe the system with the former model. Close to the critical point 
this succession is the scaling invariant hierarchy of the topology moments. This form 
of order parameter was used in [T4] to describe the glass transition in terms of critical 
dynamics. 

7. Conclusions 

The important feature of the approach developed above is the possibility of the 
sequential description of supercooled liquid. This description joins the macroscopic 
properties with the short-range and middle-range order of molecular system. One can 
conclude that the presence of a local symmetry in a real molecular system leads to the 
breakdown of the SO(3)-symmetry of the ideal continuous theory of defects. As a result 
the disclination and dislocation fields become Abelian, and this essentially simplifies 
the system analysis. In the first infinitesimal order the free disclination and dislocation 
fields are independent each other. Thus this approach is equivalent to the adiabatic 
approximation. This also simplifies the system analysis but at the same time allows to 
keep the important properties of the system. One of these properties is a decoupling 
between rotational and translational relaxation processes: according to the developed 
theory a-relaxation is the relaxation in the disclination subsystem, and /3-relaxation is 
the relaxation in the dislocation subsystem. 
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8. Appendix 1 (topology) 

The formation of the bond orientational order in liquid leads to the addition of the U 
potential, fixing the local orientation order symmetry, to the Lagrangian of the isotropic 
elastic matter. The state space of this system is the set in which the free energy potential, 
U, extremes lie. It is the orbit of the U potential symmetry grope (the ground state 
symmetry grope) on the symmetry grope, G, of the 3-d isotropic space, which coincides 
with the SO (3) rotation grope. 

The uniaxial nematic has the SO (2) ground state symmetry grope. Therefore 
its single orbit on the G is the sphere with the identified antipode points, 
RP 2 =SO(3)/SO(2) xZ 2 . The ground state symmetry grope of icosahedral, Y, has the 
three orbits on the G, which are correspond to the discrete axial rotations around the 
2-th, 3-th and 5-th order axes: SO^/Z^ xZ 2 {v\ — %v% — 3, u 3 = 5) [T5] . 

The initiation of a linear singularity (defect) in the system implies that from the 
R 3 coordinate space the R line is cut out. Thus the linear singularity presence leads 
to mapping of circle, R 3 \R 1 ~S 1 , in the M n state space. If the class of the S 1 circle 
homotopic mapping to the M n state space can not be contracted to point, ni(M n ) ^1, 
then the linear defect in the order parameter field is stable. In both cases, which are 
of interest to us the linear singularities are topologically stable, since ^(RP 2 )^! and 
7Ti(SO(3)/Y)t^1. In the non-Abelian case, when a symmetry breakdown is absent, there 
are no topologically stable defects, since 7Ti(SO(3)/SO(3))=0. 



9. Appendix 2 (the topology moment determination) 

In order to determine the disclination topology moment let us consider the tensor 
potential of an arbitrary system of the closed linear disclinations, included in the volume, 
V: 

K = e^-e\ k J pdV, (4) 

where x is the distance from the kdV disclination section to the observation point. Let 
us pick out an arbitrary point in the considered system («center of currents»), and 
signify by the R the vector from this point to the disclination section (R = x + r), then, 
provided that r -C x, in the expansion Q one can content oneself with two first terms: 



■ v R k f v ■ f r k 

,3u 4 f R k R s r s , Jjr 



(5) 



In the order of values the relation of the last pair of terms to the first one is equal to 
l/x, where I is the linear dimension of V. Therefore at x > I, when this expression is 
correct, the last two terms are smaller than the first one. Thus one can approximately 
write that 

4 = e*~4|j / IdV. 



A description of a- and ^-relaxation supercooled in liquids 



10 



However, in the case of the system of closed and screening each other topology defects, 
when 

v J kdV = 0, 

the first and the second terms in (jnj) vanish, as a result 

A\ = na 2 u (5 l a - 3e l e a ) ± + . . . . 

In this case the topology moment, Q l a = na?v (5 l a — 3e*e a ), is the basic parameter, which 
characterizes the closed linear defect system. At a great distance the tensor potential 
of this system is 

At _ Qa 

where R is the distance from the observation point to the disclination system, 
characterized by the topology moment Q . Similarly one can express the tensor 
dislocation field potential by means of dislocation topology moment, Q T . 
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